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^  Of  major  importance  in  the  determination  of  residual  stress  via  d^^ 
fraction  is  the  accuracy  of  the  measurement  of  the  scattering  angle  (2ep) 
of  a  Bragg  peak.  Ihis  determines  the  euxuracy  of  the  interplanar  (d) 
spacing  and  hence  the  strain  and  stress.  .  In  the  U.S.,  the  most  commonly 
accQ>ted  method  of  determining  peak  position  is  a  parabolic  fit  near  the 
top  of  a  peak.  (While  a  diffraction  pmk  is  not  parabolic,  this  is  a 
satisfactorjj^function  neair  the  maximui^^  The  error  in  this  procedure  has 
been  derived^Ttnd  tested,  and  it  has  been  i^wn  that  a  multipoint  fit  with 
a  least  7  points  is  rc^id  and  as  precise  or  more  prec^e  than  the  c^troid, 
the  bisector  of  the  half  width,  or  cross  correlation*^,  except  for  sharp 
peaks  in  which  case  the  centroid  or  cross  correlation  are  slightly  better. 
Thus  a  parabolic  fit  is  generally  useful  and,  since  a  least-squares  fit  to 
this  function  is  readily  carried  out  on  modem  .aicro-processors,  automation 
of  a  stress  measurement  is  possible,  including  evaluation  of  error  si  In 
this  procediire,  with  intensities  across  a  peak  at  i  2B^values,  the  var¬ 
iance  in  peak  position,  (2en)»  is: 


(2ep)  = 


(Iff)* 


iMIj). 


(1) 


2 

The  first  term  on  the  ri^t-hand  side  is  then  evaluated  for  a  parabola,  * 
and  the  second  has  been  shown  by  Wilson^  to  be: 

o*(Ii)  =  for  s  fixed  time,  t,  at  each  position,  (2a) 

»  I?/j,  for  fixed  count,  c,  at  eeKdi  position.  (2b) 


That  these  equations  are  correct  is  indicated  in  Fig.  1.  Typical  examples 
of  the  effect  of  error  in  26-  on  stress  determinations  are  illustrated  in 
Figs.  2  and  3. 

These  figures  illustrate  one  particular  goal  of  this  pamper.  There  is 
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increasing  interest  in  meeisuring  the  stresses  in  seocxid  phases  as  well  as 
the  matrix  in  multiphase  materials,  or  in  the  strengthening  phases  in  com¬ 
posites.  Because  of  their  low  volxime  fraction,  and  their  structure  and  its 
perfection,  peaks  with  reasonable  intensities  are  available  only  at  inter¬ 
mediate  angles,  where  the  errors  may  be  too  large  to  obtain  reliable  stress 
values.  Perhaps  greater  precision  can  be  Obtedned  by  curve  fitting  the 
entire  diffraction  peak  to  some  suitable  function  other  than  a  parabola. 
Furthermore,  the  speed  of  stress  measurem^ts  has  beai  greatly  reduced  by 
the  use  of  position  sensitive  detectors^'^,  by  means  of  which  an  entire 
peak  profile  is  recorded  at  the  same  time  and  no  detector  motion  is  re¬ 
quired  -  a  kind  of  digitized  return  to  filml  Why  not  use  a  fit  to  the 
entire  peak?  Ihe  data  is  already  available  in  the  same  time  it  takes  to 
obtain  the  information  for  a  peurabolic  fit. 

It  is  with  these  two  points  in  mind  that  we  have  exsunined  the  use  of 
profile  fits  in  the  me^uBurement  of  a  peak  position. 

FtKCnCNS 


Two  functions  were  chosen  to  compare  to  the  parabolic  fit.  ^e  first 
of  these  is  a  Mcxiified  Lorentzian,  proposed  by  Mignot  and  Bondot': 


I  (20j)  =  Iq[cos* 


(20.)  -  20  -6)  n 
*  P  1 


K» 


+  (2Oj-20pT  * 


(3) 


Measurements  are  made  at  many  20{  and  a  solution  is  sou^t  for  a ,  2^,  Ip 
the  maximum  intensity,  K  which  is  related  to  the  peak  width,  6  to  amount 
for  small  peak  asymmetry  and  n.  Ihe  cosine  term  forces  the  func±icai  to 
fall  more  rapidly  with  20  in  the  tail  of  the  peak  than  a  pure  lorentzian 
funcrtion.  Ihis  equation  c:an  be  mcxiified  to  include  a  K  -  K  doublet 
(with  separations): 


itotal  ^  /«i(20j)  +  %  lS(2Bi  -  A) 


(4) 


Ihe  term  6  was  chc^een  as  zero,  bec:ause  preliminary  tests  (multiple 
scans  of  peaks)  indicated  that  it  led  to  large  variations  in  20-  without 
much  improvement  in  the  fit  to  the  entire  shape.  Both  singlet  and  dcxjblet 
forms  were  tried. 

Ihe  seocxid  function  was  a  Pearscxi  lype  VII  distribution^. 


(20.  -  20„)*  „ 
■  'oC  *  — ’ 


(5) 


Ihe  terms  I^  and  a  are  scxi^t  in  the  profile  fitting  pr(x»lure.  Ihe 
parameter  m  varies  with  peak  shape  For  narrow  peaks  m  »  1  is  appropriate, 
in  which  case  Bq.  5  is  a  lorentzian.  When  m  approaches  infinity  it  can  be 
shown  that  Eq.  5  approaches  a  Gaussian  funcrticxi.  values  of  m  «  1-3  and 
infinity  were  tried.  Again,  two  terms  can  be  added  to  form  a  doublet,  and 
this  was  attended  for  m  *  3. 

Backgrcxnd  was  subtracrted.  For  the  c:ase  of  the  parabolic  fit,  this  was 
measured  at  20  (before)  the  peak  and  this  cxxistant  value  was  sub¬ 
tracted  frcxn  all  data  points.  Ihis  was  adecjuate  bec:au8e  cxily  points  in  the 

O- 


vicinity  of  the  peak  (the  tcp  15  pet)  were  employed  in  the  fit.  For  the 
other  two  fvmcticns  a  linear  variation  was  assumed  and  the  slope  and  inter- 
c^>t  were  sou^t  in  the  solution*  or  the  values  were  fixed  from  the  data. 

If  these  were  included  in  the  solution*  there  were  7  variables  for  the 
Modified  Lor^tzian  (these  two*  plus  1^*  K*  a*  n  and  2ep).  With  the 
Pearson  Type  VII  there  were  four,  (two  for  badeground  plus  1^,  and  a). 

While  a  least-equares  solution  is  possible  for  a  pardtolic  fit,  this  is 
not  ^  simple  for  the  other  two  functions  (exc^t  for  the  simpler  Pearson 
lype  VII  forms).  A  modified  sin^lex  method  was  adapted  in  this  study^'^^. 
Oonvergence  was  tested  against  the  significance  ratio  between  the  best  and 
worst  points  of  the  simplex*  whose  points  were  defined  in  terms  of  a 
goodness  of  fit*  or  reliability  index*  R:  . 

R  s  — - - 1 -  xioo.  (6) 

Z 

i=1  ‘ 

In  the  simplex  procedure,  a  multidimensional  R  space  is  fcxrmed  (n  +  1 
dimensions,  where  n  is  the  number  of  unknowns),  one  point  from  initial 
guesses  at  parameters  and  the  oth^s  from  fractional  changes  in  each  of  the 
"guesses".  Changes  are  then  made  in  these  vedues  in  a  systematic  way  to 
reduce  the  range  of  these  R  values  until  lay  some  test  it  is  found  that  all 
values  are  ess«itiedly  identical.  This  was  judged  hy  forming  r,  tlie  ratio 
of  the  reli2bility  index  of  the  wcxrst  point  in  the  simplex,  to  that  of  the 
best.  This  ratio  was  subjected  to  an  F  test^^,  that  is  r  was  calculated 
aat 


r  =  [-2-  F  + 

p,n-p,n  ‘•n-p  ^p.n-p,o 


1] 


% 


(7) 


Here  p  is  the  nuinber  of  parameters*  and  a  is  the  probability  of  incor¬ 
rectly  ooncluding  that  the  best  and  worst  points  in  the  sioplex  are  dif¬ 
ferent.  Iterations  were  continued  until  r  <  tp  „  -  .  The  test  value  was 
1.00018,  which  corresponds  to  ana  value  less^chSh^OOS.  Values  even 
closer  to  unity  were  attenpted*  but  beyond  this  value  the  error  in  2€L  was 
not  substantially  inpresved.  It  is  inportant  in  using  the  simplex  proc^ure 
to  make  appropriate  first  guesses  at  the  change  in  verifies,  so  that  the  R 
values  are  feu:  apart,  and  to  accelerate  ocnvergence  only  chafes  in  the 
parameters  which  decreased  R  were  acx^^ted;  if  an  improvement  ocoured*  the 
change  was  increased  in  the  same  direc±i(xi.  In  effect*  the  point  with  the 
worst  R  is  moved  through  the  centroid  of  the  (n  -f  1)  sided  polygon  of  R 
valies  to  lower  R.  Finally*  the  oitire  simplex  was  ocntracted  by  moving 
edl  points  in  any  iteraticxi  half  the  length  of  a  side  of  the  polygon 
toward  the  best  point.  Even  with  adl  these  precautions*  typically  nearly  5 
minutes  were  required  on  PDF  11/34  to  accxeiplish  the  ~150  iteraticxis  in 
the  cese  of  a  Modified  Lorentzian  with  7  parameters.  This  was  reduced  to 
~1.5  minutes  for  5  parameters,  and  was  "45  seconds  for  the  Pearson  Type  VII 
(with  5  parameters). 


SEKnsnao*  MiNusis 


Qjr  procedure  was  to  remeasure  a  peak  many  times  (typic:ally  10}  and  to 
examine  the  variaticxi  in  fitting  parameters.  Aocxsrdingly*  we  eqplcsyed 
various  statistical  tests  to  jui^  the  results.  The  mean  of  2ep  and  R  for 
any  set  (n)  of  peaks  fit  \xy  a  single  method  was  obtained,  as  well  as  the 
ebserved  variance*  S*.  ^ 


(8) 


S*  (<2®p>)  =  nTT  f  -  <28p>)*  . 

cuid  as  well  the  true  variance  »  S^/n.  Ihe  range  around  the  mean  with 
95%  confidence  limits  was  formed  as  <2ep>  +  1.96a  (2ep).  Ihe  oon- 
fidance  limits  for  the  variance  was  established  with  a  chi-squared  test 
(C  )  at  the  1-0  confidence  level  as  follows: 


o* 


<o 


* 


=  (9) 


fbr  a  given  peak,  the  use  of  different  profile  fitting  methods  could 
lead  to  different  2^.  A  test  was  made  to  ascertain  if  such  differences 
were  read,  or  due  to  hunting  statistics.  Let  i  be  the  number  of  different 
profile  methods  for  any  one  peak.  Ihen  the  “profile  method"  variance  was 
defined  as: 


S',  «V' '  *- 


2  (2e.„ 

i 


-  <2V>’ 


(10) 


where  <2ep>  =  j  ^  number  of  peaOc  positions  for  the 

different  methods  for  the  same  peak. 

The  pooled  variance  for  adl  methods  is  defined  as: 


(11) 


With  n  peaks  for  each  method,  there  are  ttven  (n-1)  degrees  of  freedom 
for  each  method,  and  ^  (n-1)  for  (and  -1]  for  S^.  The  ratio: 


(12) 


was  then  tested.  If  Bq.  12  is  less  than  the  tabulated  F  value  for  some 
confidence  level,  the  difference  in  mean  position  is  due  to  random  counting 
statistics  (with  that  confidence). 
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Sandies  were  chosen  to  represent  a  topical  range  of  peak  intensities 
and  peak  widths,  and  to  include  second-phase  peaks.  Their  preparation  and 
the  peak  characteristics  are  described  in  Table  I.  The  diffraction  cx»di- 
tions  are  in  T^le  II.  The  data  were  obtained  on  a  micrc3proc»ssor  ocxitrol- 
led  diffractometer.  The  parabolic  fit  was  performed  on  line,  and  data  were 
chtained  in  a  three  stage  point-oountii:g  prcx^iure  described  in  Ref.  2,  to 
IcTcate  the  peak  and  the  range  of  26  cxsvering  the  tc3p  15  pet.  Seven 
points  in  this  range  were  employed.  For  the  other  functicxis,  data  was 
obtained  with  the  same  system  (without  removing  the  specimen)  at  .02926 
intervecls  (each  cxxmted  for  15  secxsnds)  across  the  entire  profile,  and 
transferred  to  a  FDP  11/34  minicximputer  for  data  processing.  The  data  were 
not  ccrrected  for  the  Lorentz-polarizatim  fac±or  or  scattering  factor 
variaticxi,  as  wcxild  be  needed  in  ac:tual  stress  analysis;  only  the  fit  to 
the  shape  and  the  value  (and  error)  of  26  were  of  interest  here. 
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.  Ab  with  the  parabolic  fit,  the  Pearson  Type  VII  distribution  fits  only 
the  main  part  of  the  peak,  but  falls  more  rapidly  in  the  tails  than  the 
observations,  the  difference  being  larger  as  "m"  increases.  (The  region 
fit  was  actucilly  similar  to  that  for  a  parabola.)  On  the  other  hand,  the 
Modified  lorentzian  fits  the  entire  shape  quite  well.  A  comparison  of  the 
various  fits  with  ten  recordings  of  three  different  peaks  is  given  in  T^le 
III,  and  samples  of  the  fits  are  shown  in  Fig.  4.  A  fixed  lineeu:  back¬ 
ground  was  employed  for  the  results  in  this  table;  solutions  with  a  var¬ 
iable  background  had  worse  errors  in  peak  position,  although  with  the 
•significance  tests  at  the  95  pot  level  these  e^^parent  differences  were  not 
necessarily  real.  More  importantly,  the  errors  in  the  fitting  pareuneters 
showed  drastic  decreases,  and  the  values  of  K  and  I^  epproached  the 
measured  peak  width  and  intensity  when  the  background  was  fixed  prior  to 
the  solution  for  the  other  parameters. 

It  is  evident  from  the  table,  that  the  Modified  lorentzian  and  peura- 
bolic  fit  provides  similar  error  vadues  for  £dl  peak  shapes.  This  is  true 
of  the  Pearson  Type  VII  function  as  well  but  the  value  of  m  is  different 
for  each  peak  type.  Thus,  some  knowledge  of  the  peak  shape  would  be 
required  prior  to  fitting.  Which  would  make  automation  more  difficult  than 
for  the  other  functions.  Also,  the  Goodness-of-Fit  values  are  quite  high 
with  this  function. 

Comparing  the  method  and  pooled  variances  indicated  that  the  sliest 
differences  in  peak  positions  in  Ihble  III  for  the  various  functions  are 
significant  in  all  cases  except  for  the  broad  peaks.  In  this  latter  case 
the  low  intensity  leads  to  larger  scatter  so  that  any  difference  is  masked. 

Note  especially  that  the  parabolic  fit  gives  the  lowest  errors  of  any 
method  for  weak  or  broad  peaks. 

The  60/40  breuss  sanple  exhibited  a  doublet  whose  resolution  varied  with 
tilt  of  the  san^le  to  the  x-ray  beam  (as  would  be  done  in  a  stress  measure¬ 
ment).  Exanples  of  fits  to  this  peak  are  given  in  Fig.  5,  and  a  summary  of 
results  in  Thble  IV.  Ihe  singlet  fcxrms  of  the  Modified  Lorentzian  and 
Pearson  Type  VII  functions  gave  values  of  A  20C((J»*  0°)  -  ((|»»  45°)] 
closer  to  the  valxies  of  the  parabola,  despite  the  fact  that  the  <R>  values 
were  lower  with  the  doublet  form.  From  the  analysis  of  variance  at  both 
angles,  the  differences  in  peak  positions  with  each  fit  technique  in  this 
table  were  significant. 

Possible  limitations  in  the  data  that  could  ocxmr  in  prac:tic:e  were  also 
e;q>lored.  Firstly,  it  may  not  be  possible  to  record  the  entire  peak; 
another  one  nearby  may  overlap  on  one  side,  or  the  equipment  itself  could 
preclude  recording  the  entire  peak.  Some  results  are  shown  in  Table  V. 

With  a  Modified  Lorentzian  and  a  sharp  peak  the  "correct"  (parabolic)  peak 
position,  a  low  error  in  this  value,  and  a  low  R  value  are  all  obtained 
with  data  that  only  just  reaches  the  peak's  maximum.  Simileu:  tests  with  a 
broad  profile  showed  that  a  wider  range  was  necessary,  but  only  a  few 
tenths  of  a  degree  2e  beyond  the  maximum  is  adequate.  Analysis  of 
variance  tests  confirmed  that  beyond  68.82°,  any  change  in  peak  position  is 
solely  dxie  to  random  fluctuations,  and  the  peak  values  obtained  with  data 
ip  to  and  beyond  this  angle  are  not  distinguishable  from  the  value  for  the 
entire  peak,  at  95  pet  cxxifidence. 

Next,  with  the  Modified  Lorentzian,  the  effect  of  time  per  data  point 
and  the  nunber  of  data  points  were  explored,  Tbble  VI.  The  error  is  not 
statistically  different  at  the  shorter  counting  time,  or  smzd-ler  number  of 


[For  further  details  the  reaier  is  referred  to  Ref.  12.3 
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'^'Ihe  parabolic  fit  has  the  best  overall  ability  to  determine  peak 
positions  over  a  wide  range  of  shapes. 


A  Modified  Lorentzian  can  be  quite  helpful  if  only  a  part  of  a  peak  can 
be  esqslored,  especially  for  sharp  peaks. 


Ihe  errors  in  a  Modified  Lorentzian  fit  are  not  very  sensitive^ to 
counting  time.  With  a  positi<»  sensitive  ^detector  precision  comparable 
to  that  for  a  parabolic  fit  can  be  obtain^  in  about  one  tenth  the 
time.  As  adl  of  the  data  are  a^ecorded  for  either  fit  with  this  type  of 
detector,  a  considerable  saving  in  measurement  time  is  possible  with  a 
Modified  Lorentzian  function.  Unfortunately,  the  fit  itself  takes  the 
order  ^  1  minute  on  a  miniprocessor.  If  this  time  could  be  reduced 
with^ specially  designed  microprocessor  for  this  purpose,  this  fun- 
ctjion  offers  a  way  of  dreistically  reducing  the  time  for  stress  measure- 
its  in  the  field,  beyond  that  alreai^  adnieved  by  the  use  of  a  PSD. 


ss-of-Fit  is  not  necessarily  an  accurate  gage  of  error  in  peak 
positioiw  (Ebr  exaople,  in  Thble  III,  part  (c),  similar  values  of  R 
are  assmnted  with  widely  different  peak  locations.) 
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TSBEZ  II 

OPERATING  GCNDITICNS 


Divergent  Slit 


Receiving  Slit 


Annealed  1008 
1074,  Brass 


Beam  Size  cn 
at  =  oP 
(approximate) 


2.^nm  X  2.5inm 


Tuibe  Target 


Tube  Voltage- 
Current 


40kv 

35kv 

40kv 

35kv 


Annealed  1008 
Deformed  1008 
1074 
Brass 


Filter 


Vanadium  Oxide 


Soller  Slits 


None  used 


uans. 

QUITS 


TtW  (SeCONOSl 


FIG:  1.  Comparison  of  statistical  error  in  intensity  measurements  with 
observed  error.  Left:  Fixed  counts,  Eq.  2a,  Right:  Fixed  time 
measurements,  Eq.  2b.  Twenty  five  replications  of  a  point  on  a 
211  ferrite  reflection  from  a  1008  steel  was  employed  with 
filtered  CuK^  radiation. 


nABlE  III 


OMPARISON  CF  FIT  MEJIHQDS^ 


MEHIiOD 

MEAN 

^  P 

OBSERVED 

ERROR 

(xl0"2 

degrees) 

95%  LIMITS 

ON  ERTOR^ 

(xl0"2 

degrees) 

MEAN 

QOMMESS-OF- 

FIT 

(Annealed) 

Fexrite  110, 

1008  Steeli 

68.82193 

.287 

.197,  .514 

29.91 

Pearson  TVpe  VII^ 

m  =  3 

68.81368 

.473 

.326,  .863 

30.24 

m  -  2 

68.82818 

.172 

.118,  .314 

32.43 

m  -  1 

68.82902 

.282 

.194,  .515 

41.77 

m  s 

68.80705 

.172 

.118,  .314 

33.43 

Modified  Lorentzian^ 

68.83298 

.135 

.093,  .246 

1.75 

B)  Curve  Type:  Broad 

(deformed) 

Fbrrite  110, 

1008  Steel 

Pcurabolic 

68.53719 

.173 

.199,  .316 

1.16 

Pearson  Type  VII 

m  =  3 

68.53853 

.511 

.352,  .870 

20.53 

m  s  2 

68.54362 

.334 

.229,  .609 

9.95 

m  *  1 

68.53S95 

1.40 

.963,2.56 

17.07 

m  * 

68.54447 

.277 

.191,  .506 

11.21 

Modified  Lorentzian 

68.54165 

.351 

.242,  .641 

1.88 

C)  Curve  Type:  Weak 

Pe^C  112  + 

021,  1074  Steel 

Parabolic 

57.42794 

1.05 

.723,  1.92 

2.79 

Pearson  TVpe  VII 

m  =  3 

57.35688 

.263 

.181,  .480 

26.5 

m  =  2 

57.36416 

1.78 

1.22,  3.25 

17.2 

m  =  1 

57.36556 

.785 

.540,  1.43 

6.98 

m  = 

57.40937 

5.15 

3.54,  9.40 

4.49 

Modified  lorentzian 

57.40813 

1.28 

.881,  2.34 

2.11 

TI  15  Seconds  per  point,  10  peaks  per  mean.  Fixed  background  -  5  points 
either  side  of  peak  used  to  fix  line  (for  parabola,  see  ref.  2). 

2.  Parabolic  fit  is  7  pt.  fit  to  top  15%.  A  sii^le-valued  background  is 
subtracted  before  fit. 

3.  All  Pe^u:son  I^pe  VII  and  Modified  Lorentzian  utilized  90  pts.  taken 
at  .02°  2  increments,  except  for  FejC,  for  v<Aiich  58  points  were 
measured. 

4.  Values  listed  are  low  and  high  limits  respectively. 


60/40  Brass,  211  peak,  Crl^  Radiation,  10  sccin  replications,  78  points  per  pecdc). 


TABLE  V. 


PARTIAL  FIT  OF  SHARP  PEAK  (1008  STEEL)  TO  MODIFIED  LOREMTZIAN 


TWO  THETA  <20_>  (BSERVED  <R> 

VALUE  OF  ,05!  X  ERROR 

lAST  POINT  ' 


68.70 

68.77712 

.296E-1 

1.48 

68.74 

68.79541 

.5a4E-2 

1.47 

68.78 

68.82561 

.341E-1 

1.70 

68.82 

68.83322 

.lOlE-2 

1.74 

68.86 

68.83303 

.836E-3 

1.74 

68.90 

68.83298 

.135E>-2 

1.75 

68.94 

68.83226 

.113E-2 

1.75 

68.98 

68.83322 

.llOE-2 

1.74 

69.02 

68.93300 

.937E-3 

1.72 

iO  second  count  per  point,  0.02®  2e  intervcLLs. 

O 

Lowest  two  theta  veLLue  in  all  cases  was  67.50. 

'Value  of  Last  Point'  indicates  the  highest  two  theta  veilue  used  in 
analysis. 

Peak  value,  ais  determined  from  the  full  angiilar  range  fit,  is  68.83299 
±  .00135  “20. 

A  fixed  background,  obtained  from  10  data  points  on  the  low  angle  side  of 
the  curve,  was  used. 


TARTJt  VI 

EFFECT  OF  NUMBER  OF  POIOTS-MODIFIED  LOREOTZIAN,  BROAD  CURVE  (1008  STEEL) 


TIME 

PER  POINT 
(SEC.) 

NUMBER  OF 

POINTS^ 

MEAN  PEAK 
POSITIC*! 

(°  20) 

OBSERVED 

ERROR 

(xl0“2) 

<R> 

1 

22 

68.54718 

.572 

3.09 

90 

68.54763 

.434 

3.59 

15 

22 

68.54224 

.360 

2.07 

90 

68.54165 

.351 

1.88 

1.  The  22  points 

were  selected  as 

every  fourth  point  of  the  full  90  point 

fit.  Background  line  determined  from  5  points  on  either  side  of  the 


Estimated  statistical  error  in:  a)  "d"- 
spacing  and  b)  stress,  as  a  function  of 
diffraction  angle.  A  seven  point  para¬ 
bola  was  assumed. 
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